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NOTATIONS 



t — time 

e — the basis of natural logarithms 
j = /- 1 

Hi - the mass of the Earth 
m' — the mass of the Moon 
m" — the mass of the Sun 

k' , /j.' , v‘ — the equatorial components of the unit vector directed from the 
center of the Earth toward the Moon 

r' — distance of the Moon from the center of the Earth 

a' — mean value of r' defined in such a way that the constant term in 

the trigonometrical expansion of a'/r' is equal to 1 

k", fj." , v" - the equatorial components of the unit vectors directed from the 

center of the Earth toward the Sun 

r" — distance of the Sun from the center of the Earth 


a" — mean value of r" defined in such a way that the constant term in 
the trigonometrical expansion of a"/r" is equal to 1 
a — semimajor axis of satellite's orbit 

i - the inclination of satellite's orbital plane toward the equatorial 
plane 

Q - the right ascension of ascending node of satellite's orbital plane 
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Qj - the mean motion of the ascending node of satellite's orbital plane 
e - the eccentricity of satellite's orbit 
n — the mean motion of the satellite 
oi — the argument of the perigee of satellite's orbit 

7T — CO + fi 

M — the mean anomaly of the satellite 
R - the equatorial radius of the Earth 
a = R/a 

a' = R/a' — the lunar parallactic factor 
a" = R/a" - the solar parallactic factor 
-t, I', F, D, T — the arguments of the lunar theory 

, l‘ v Fj , Dj , Tj - the mean motions of the lunar arguments 
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ON THE TIDAL EFFECTS IN THE MOTION OF 


EARTH SATELLITES AND THE 
LOVE PARAMETERS OF THE EARTH 


INTRODUCTION 

At the present time, the tidal effects in the motion of artificial satellites 
attract considerable attention from specialists in geophysics and celestial 
mechanics, mainly because from these effects one can determine the elastic 
properties of the Earth as they are observed from extraterrestrial space. 

The disturbing potential pertaining to the problem is obtained as the analytical 
continuation of the tidal potential from the surface of the earth into outer space. 
This continuation must include parameters which characterize the earth’s elastic 
response to tidal attraction by the moon and the Sun. There are in use at the 
present time several versions of the exterior tidal potential. 

In setting a priori the tidal potential in the form of a functional linear in 
the elastic impulse responses of the earth we are being guided by our previous 
expansions and also by earlier works on the theory of tides. In our previous 
work (Musen and Estes, 1971), (Musen and Felsentreger, 1972) we expanded 
the exterior tidal potential into a sum of products of spherical harmonics in 
Maxwellian form. The first factor in each product depends upon the rectangular 
coordinates of the satellite and the second upon the coordinates of the disturbing 
body. 
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The elastic parameters of the earth were introduced in the form of factors 
(Love numbers) attached to harmonics in the expansion of the potential. The 
simplest assumption about the Love numbers is that they are the same for all 
spherical harmonics of a given degree. Such an approach was taken by Kozai (1965), 

Fisher and Felsentreger (1966), and Smith, Kolenkiewicz and Dunn (1971). As 
we stated before (Musen and Felsentreger, 1972) this approach is equivalent to 
introducting "global" Love numbers obtained by a double averaging process, as 
applied to the local Love numbers, performed over the whole earth. 

A different approach was suggested by Kaula (1969). He pointed out that 
the dependence of the local Love numbers on latitude can produce noticeable 
long period effects in the motion of a satellite, whereas the dependence on longi- 
tude will produce predominantly small short period effects of the period of one day 
or less. If we are interested only in the long period tidal effects in the motion 
of a satellite, then the dependence of Love parameters upon latitude shall be 
taken into consideration before we start the analytical continuation of the tid^l 
potential into outer space. Neglecting the dependence upon longitude is equivalent 
to averaging the local Love parameters along the parallels, and not ove 1 ' the 
whole earth. The resulting expansion of the exterior tidal potential is, as before, 
a sum of products of spherical harmonics. Each pair of harmonics of a given 
degree and order has the identical Love number. 
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It is of interest to note that this idea can be found already in the theory of 
tides by Laplace. Among the new approaches the one developed by Munk and 
Cartwright (1966) deserves special attention. The tides are represented by a 
sum of convolutions between the spherical harmonics and the earth's impulse 
responses. The impulse responses replace the Love numbers. The earth's 
tide is a very complicated process. Continuous masses are moving and the 
mutual influence of adjacent (in time and space) configurations takes place. The 
evidence of such "cross-effects" is supported by the observational material. 

A theory should be devised to include the influence of these combined cross- 
effects on the motion of an artificial satellite. It is clear that the resulting 
potential cannot be obtained from the "first principles" of celestial mechanics. 

There is a time lag in the exterior tidal potential acting on the satellite. 

This fact, together with the possibility of interpreting Love numbers as Dirac 
impulse functions leads to a representation of the tidal potential in the form of 
a convolution between the Love impulse responses and the spherical harmonics 

depending upon the coordinates of the moon (Sun). This is a general form of the 
exterior tidal potential which includes the influence of past configurations on 
the motion of the artificial satellite. The impulse responses which replace the 
delta functions mentioned previously are, in all probability, continuous, fast 
decreasing functions. They are determined from the condition that the 
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mean quadratic error between the observed and the computed tidal effects 
shall be a minimum, which leads to a system of Fredholm integral equations. 

The transformation of these equations from the time domain to frequency domain 
provides us with simple relations which can serve as a check for the existence 
of Love numbers. 

We formulate the integral equations and the corresponding relations in the 
frequency domain for the main problem only, considering harmonics of the second 
degree in the tidal potential. In the present exposition we resort, as in our 
earlier work, to an expansion into trigonometric series with arguments 
V , F, D, F of the lunar theory and Q , n of the satellite. 

TIDAL DISTURBING FUNCTION 

Let A.' , fj.' , v' and\" , /j." , v" be the equatorial components of the unit 
vectors directed from the center of the Earth toward the Moon and the Sun, 
respectively. 

Let r' and r" be their distances from the center of the earth and a' f a" 
be the mean values of r' and of r" , respectively, defined in such a way that 
the constant terms in the trigonometrical expansions of the parallaxes a' /r' 
and a"/ r" are equal to 1. 

The basic spherical functions which appear in the expansion of the exterior 
lunar and solar tidal disturbing functions associated with the "main problem" 
are (in Maxwellian form): 
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for the "lunar" disturbing function: 


c o ( t )=( 7 7] (1-3A 


C +1 Ct) = (77) Cm' + i *•'). 


C -l C*) =(^7 Cm' - j *') 


(!’) 


c + 2 (t) =(77) CM' + j V) 2 


c -2 Ct) =(77) (M' - j ^') 2 


and, similarly, for the "solar" disturbing function: 


c o Ct) = — (1 -3/ 2 ), 


<1 (t) =(~) Cm" + J *'), 


c — 1 (t) = (- ) Cm" - j V'), 

1 1 * 


<2 ^> = (777) Cm" + j A.") 2 , 


(l") 



Cm" - j V') 2 
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Introducing the elastic impulse responses of the earth, 

W 0 (*)• W 1 W 2 (*)• 

we can rewrite the exterior tidal disturbing functions acting on the artificial 
satellite as given in our previous work (Musen and Estes, 1971) in the following 
compact form: 


V' = n 2 a 2 k' |^1 - ^ sin 2 ij (c^ * w Q ) 

+ [e + i^(c^ 1 *w 1 )+e _i ^(c^ 1 *w 1 )]sinicosi (2 1 ) 

- [e +2 j 0 (c^ 2 * w 2 ) + e -2 > ^ (c '_ 2 * w 2 )] s in 2 i| (1 - e 2 )“ 3/2 , 

V" = n 2 a 2 k" ^1 - sin 2 i \ (c" * w Q ) 


where 


+ [e + j ^ (c" x * Wj) + « > ^ (c"j * Wj)] s in i cos i 
- [e + 2 ' ^ (c " 2 * w 2 ) + •©” 2 > ^ (c " 2 * w 2 )] s i n 2 i\ (1 - e 2 ) _3/2 (2") 


m 

IT 


of a 


1 3 


k" =— a? a " 3 

ftl 


and the asterisk designates the convolution, taking only the past into account: 


* 



c 


n 

k 



(t - T) WL 


(t - T) W[ k 


(t) d r, 


(r) d r 


k = 0, ±1, ±2. 
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For combined luni-solar tidal effects the disturbing function becomes: 


V = n 2 a 2 k‘ 


1 - o sin M ( c o * V 


+ [e + i ^ (c +1 * Wj) + e i ^ (c_j * w x )] s in i cos i 

_ [« +2 i 0 ( C+2 * w 2 ) + o“ 2 (c_ 2 * w 2 )] sin 2 ij- (1 - e 2 )" 3/2 , (2) 

where we set 


and 


C - C 1 + K c" 
s s s 


K 


m 




All the functions 


c ' . c" s , C s (s = 0, ± 1, ± 2) 

are expanded into Fourier series with numerical coefficients and arguments 

Z, V , F, D, r, n and n . 


The details of this expansion (by electronic computer) are described in our 
previous article already quoted. In (2’) - (2) we retained only the long period 
and the secular terms. The short period terms, with periods equal to the period 
of revolution of the satellite or less, are removed by means of averaging the 
disturbing function over the instantaneous orbit of the satellite. In fact, we 
are interested only in the long period tidal effects in the elements of the 
satellite and therefore eliminate the secular effects in n and fi by omitting 
the constant terms in the Fourier expansions of , c" and c Q . 
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We shall continue, however, to use the same notations, c' 0 , c", c Q , for the 
corresponding series with the constant parts omitted. We set: 


a' (t) = J « + i sQ c' (t) d t = a' (t) + j yS; (i), 
a" (t) e + i sfi c" (t) d t = a" (t) + j /3" (t). 


(3) 


where 


a s (t) = J c s (t) d t = a s (t) + j /? s (t). 


a s = a s + K a s’ + K 

0. = - z 3 


a = a , 

S -s’ 


From (l 1 ) we deduce: 


a +i _ a -i 


r ' a '\ 3 

— -j v' (/x' cos Q. - sinfl) d t, 


+ 2 



7 ) [(At' 2 - \' 2 ) cos2fi-2^>' sin2 0] d t, 


iL) [(^'2 _ x.'2^ sin 2 n + 2 /x' cos 2 fl] d t. 


(4) 


8 



Similar expressions can be deduced for a" x , /3^, a " 2 and . The 
integration of Fourier series in (3) and (4) is performed in a purely formal 
manner, without adding constants of integration. As a result, Fourier series 
for a^, a", a s contain only purely periodic terms. 

We deduce after some easy transformations: 

(c' * w|, s |.) d t = a' * (w| s |,e jsfilt ) , 

Je jsfi (c"*w| s |)dt = a;'*( w|s| .e jsn ^), 

(c s * w )# |.) d t = a, * (wj ,| e ^ 

As before the convolutions are computed taking only the past into consideration 
and the integration with respect to time is performed in a formal manner. 

PERTURBATIONS IN THE ELEMENTS DUE TO TIDES 
Making use of the Lagrange equations 

d S 1 1 d V 

d ^ n a 2 / 1 - e 2 sin i ^ U’ 


d S Q 

d t 


= + 


1 

n a 2 / 1 - e 2 sin 


dV 
d i’ 
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d S 77 ^2 3 V /l - e J 3 V 

— -j. 

d * n a 2 / 1 - e 2 ^ * n a 2 e ^ e 


d 8 M _ _ 1 - e 2 d V _ 2 / 3V 
d t n a 2 e 3 e n a 2 V 3 a 

We deduce for the tidal perturbations in the frame of the main problem: 

d ^ * = - n K j { [e + j^(c + *w 1 )-e - ^(c * Wj)] cos i 

dt (1 - e 2 ) 2 

- 2 [-e + 2 > ^ (c +2 * w 2 ) - e~ 2 i ^ (c_ 2 * w 2 )] sin i} 



d § n 

d t 


= + {- 3 (c * w ) cos i 

(1 - eV 

+ [e + (c +1 * Wj) + e _i ^ (c_ 1 * w 1 )] (cosec i - 2 sin i) 


(7) 


-2 [e + 2 i ( c+2 * + e 2 i ^ (c_ 2 * w 2 )] cos i} , 


d S 7r n k' f 3 , c 

d t (1 _ e 2 ) 2 l 2 


cos 2 i - 2 cos i - 1) (c Q * w Q ) 


+ (5 cos 2 i + 3 cos i - 1) tg - [e + ^ (c +1 * Wj) + ■& j ^ (c_j * Wj)] 


- (5 cos i + 3) (1 - cos i) [e + 2j ^ (c +2 * w 2 ) + e 2i ^ (c_ 2 * w 2 )]j- 


(8) 
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3 


d S M _ 3 

d t + ^2 


V 


n /< 


_ p 2 n 3/2 


(1 - C 2 ) 

[e + j n (< 


1 -Isin 2 i) (c Q * w Q ) 


j) + -e j ^ (c_j * Wj)] sin i cos i 


+ 1 


- [« +2j ^ (c 2 * w 2 ) + e 2j 0 (c_ 2 * w 2 )] sin 2 i|- 


After the integration we obtain: 


8 i = - 


n k j 


(1 - e 2 ) 


2x2 


- 2 


a +2 * \ W 2 6 


{[ a+ i * (wj - a., * (w, e ,fl ‘ ‘)] cos i 

(w 2 6 +2jfil *) - a_ 2 * (w 2 e 2) filt j]s in i} , 


( 10 ) 


S n = + - - {- 3 ( a o * V cos 1 

(1 - e 2 ) 2 1 

r / 

+ [a +1 e / + a_j * ( Wj « 1 (cosec i - 2 sin i) 


- 2 


[ a +2 * ( w 2 


w„e '1 +a _ 2 * (w 2 e‘ 2ifllt )] cos i} , (11) 


n K » J (2 

S 77 = + ■(— (a. * w ) (5 cos 2 i - 2 cos i - 1) 

(1-e 2 ) 2 v.2 ° 

+ [ a +1 * (wj « 1 ) + a_j * (wj e 1 ) ] (5 cos 2 i + 3 cos i - 1) tg L 

- 2 [a +2 * [w 2 e 3 1 ) + a_ 2 * (w 2 e 2jfil (5 cos i + 3) (1 - cos i)j (12) 
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S M 


n k 


{(l ~f sin2 A ) ( a o * w o' 


(1 - e 2 )3/2 

T *( +iCl A *( • • 

3 +l * \ W 1 e / + a -! * \ W 1 ^ /J Sin 1 C ° S 1 

[a +2 * (w 2 e +2)Dlt ) + a _ 2 * (w 2 « 2)filt )] sin 2 i} 


INTEGRAL EQUATIONS FOR THE IMPULSE RESPONSES 

We determine the impulse responses from the condition that the mean quadratic 
error between the observed and the computed tidal perturbations in the elements 
shall be a minimum. The application of this principle leads to a system of 
Fredholm integral equations of the first kind for the impulse responses. As we 
stated before all functions with which we operate in the computation of the tidal 
effect and in the formation of the integral equations are given in the form of 
Fourier series with numerical coefficients. 


Such a form facilitates greatly the computation of mean values of the 
form 

E[f(t)] = li m f f ( t ) d t 

T - * 00 2 T J. T 

and of the cross- and autocorrelations 

a p,q c T ) = E [a p (t + r) a* (t)] 

= E t 3 p (t) 3 q (t - T )1 > 

which we introduce in the exposition. In the process of computation of a pq (r) 
only terms with the identical arguments ina (t) and a* (t) produce a non-zero 
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term in the cross-correlation. One of us (R. E.) developed a program for the 
computation of a pq in the form of Fourier series. The arguments in the output 
are linear combinations of the arguments V , F, D, P , tt and ft, but with t 
replaced by r and with the constants of the phase omitted. For example, making 
use of (10) we deduce integral equations for w x and w 2 by minimizing the 
quadratic functional: 


J = E 

- 2 
where 


r 

b i + j s 

IT + j ft , r - j ft . fi 

{[f+i (t - t) « -a_ 1 (t-r)e ‘J 

Jo 



1 w i 


(r) cos i 


a +2 < t “ T ) 


+ 2j ft,i 


-2j fij 


- a -2 (t - T ) « ' 1 ] W 2 (r) sin ij- dr 


(14) 


(1 - e 2 ) 2 

and S i (t) are the observed tidal perturbations in inclination. After some easy 
transformations and taking 


< (t) = a - k (t) 

into account, we obtain assuming that the order of 
can be interchanged: 



/•CO 

= E | S i | 2 

+ 2 s [T 10 (r) Wj (r) cos i 



r 

/*CO 

-■* f 

[£„ (T, V) Wj (r) Wj (i/) cos 2 

-'o 

J o 


v) Wj (r) w 2 (v) sin i cos i 

+ -t 22 (T, V 

) w 2 ( r ) w 2 C 1 ') sin 2 i] d r d v 


averaging and of the integration 


+ ^20 ( T ) W 2 c T ) sin i] d r 


( 15 ) 
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Taking the gradients with respect to w x and w 2 and setting them equal to 
zero, we obtain: 


*10 (T 


/% +00 

) - S L '& 11 (T, V ) Wj (v) CO s i - l 12 (r, v) w 2 (v) s in i] dv =0 (16) 

J n 


J* <x> 

[- -& 12 (^, T) w x (v) cos i + l 22 ( t , v) w 2 (v) sin i] d v = 0, ( 17 ) 


As we will see later 

*i 2 ^ T ’ ^ S °* 

and we obtain a pair of uncoupled integral equations 

l lQ (r) -s cos i -t u (r, v) w x (v) d v = 0 

Jo 

*20 - s sin 1 * 22 ( r ’ w 2 (v)dy=0 

•Jn 


(17’) 


(18) 


(19) 


It is of interest to note the presence of the orbital inclination in these equations. 


where 


explicit form of 

*io 

(r). 

*20( T )» l ll T ’ V 

). 

* 12 ( r > v ) 

and 

*22 (r. 

v ) are: 

*io ( T ) 

= J 

t 

+ j A. T 

« k + i ( r ) ' 

- -e 

■ jf V k 

k -i 

( T >_ 


(20) 

*20 ( T > " 2 


r - 

i-2 j n.r 

_ , 

2 j n,r 

1 


(21) 

j 

It e 

1 k +2 (r) + 

e 

1 k -2 


i 


k -l 

( r ) 

= E [§ i (t) a_ t 

(t 

- T)] , 



(22) 


k +l 

( r ) 

= E [8 i (t) a +1 

(t 

- r)] , 



(23) 


k -2 

( r ) 

= E [8 i (t) a _ 2 

(t 

- r)] , 



(24) 


k+ 2 

( r ) 

- E [S i (t) a +2 

(t 

- T )] - 



(25) 
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and 


Z n (r, v) 


+j n. (t + i/) 

+ e 


a + l,-l ( T - v ) 


-j n, (r- v) 

- ■& 

+ j n, ( t - v > 

- -e 

-j n, (r + v) 

+ e 


a + 1 , + 1 ( r - v ) 
3-1,-! ( r - v ) 
a _l, + l (T - V) 



12 


( r > 


+ jfij(r+2v) 

Z') = + « 

-j n, (t- 2v) 

- e 

+j n. ( r- 2v> 

- e 

-j 0. (T+ 2v) 

+ e 


a + 2 , - 1 C T * V ) 
a + 2 , + 1 ( T ~ ") 
3 - 2 ,-! (T - V) 

3 - 2 , + 1 ( T " V > 



+2jn,(T+v) 

v) - + -6 1 

-2j H, (T-v) 
- •© 


+ 2 j n. ( r - v ) 
•e 


-2jfi l( r + v) 
+ « 1 


a 


+ 2 , 


-2 


8 + 2 , + 2 


a 


- 2 , 


-2 


3-2, +2 


(T _ V) 
(r - V) 


(T - V) 


(T - V ) 


Substituting 

\ ~ a k + j @k’ k = + 1, ±2 

into (26)-(28) and defining: 

p n (r) = E [aj (t + r) aj (t) - P x (t + r) ^ (t)] , 

q n (r) = E [a % (t + r) (t) + P x (t + r) ( t )] , 

r n (t) = E [a x (t + r) (t) + (t + r) /3 } ( t )] , 

s n (r) = E [aj (t + r) ^ (t) - (t + r) a x (t)] , 


( 26 ) 


( 27 ) 


( 28 ) 


( 29 ) 


15 




^22 ( T ’ v ) = 8 [p 2 2 ( T ~ cos ( T + 2 V) - q 22 (r - V) sin Qj (r + 2 v) 

- r 22 ( T ~ v ) cos ( T - 2 v) - s 22 (r - v,) sin fij (r - 2 ^)] (34) 

The algebraic computations performed on the machine by one of us (R. E.) 
yield 

p u ( r ) = da ( T ) = p 12 ( T ) = q 12 ( T ) = r i 2 ( T ) = s i 2 ( r ) 

= P 2 2 ( T ) = d 2 2 ( T > = r 2 2 ( T ) = 0 

and the equations (32) -(34) become 
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'fq 1 c T - v ) ~ - 2 r H ( T - v), COS Qj (r - v) + 2 s n (r - u) sin fij (r - v), 
^12 (t, v) = 0 , 

^22 ( T ’ = - 8 s 22 ( T ~ v ) sin ( T - 2 V). 


The equations (18) -(19) can be put into a simpler form by transforming them 
to the frequency domain. In the actual computations with the assumed numerical 
accuracy the kernels of the integral equations, -C (r,v ) and (r, v), also 
^ 2 0 ( r ) * ^i 0 ( r )> as wel1 as a pq ( T _v )> a p (r) (P = °» ±1 » ±2 ) are trigonometric 
polynomials with arguments linear in r and v and, consequently, they are 
almost periodic. In transforming an almost periodic function g(t) to the frequency 


domain the averaging operator E takes the place of the Fourier transform, 


G (f) = E t [e“' ft g (t)] = lim ~ f 

j — oo ~ T 


g (t) d t, 


and G(f) is different from zero only if f is a Fourier exponent of g(t). 

We are not presupposing, however, that Wj (v), w 2 (^) are almost periodic. 
In the transformation of (18)-(19) to frequency domain the standard Fourier 
transforms 

/% +C 0 

w i (f) = 2 V fv w. (v) d v 

^ ♦'-CO 

will appear in the final result. 
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From (26)-(28) it is evident that the integrals in the equations (18)-(i9) are 


sums of terms of the form: 


h m,n (t) 
p , q v ’ 



(m r + n v) 

a p.q ( T ~ V ) W O) d V 


w = Wj, w 2 . 


Taking into account that h^" (r) is an almost periodic function of r and 
assuming that the operators of averaging and of integration commute, we have: 


H m ’ n (f) 
p, q v ' 


r co 

I r -irjf-mQ,) "1 

= E rL e * 


+ i n Q, v 

«• W (v) d v 


and, after the change of variable under the E T -sign, 


< f > = 2* E r a (r)] 


a P ,q < T >J W [f - (m + n) Qj] 


f-jT(f-ma) *\ 

- 2 7T E r je E t [a p (t + r) a_ q (t)] W [f - (m + n)fi x ]j- 

where W(f) is the Fourier transform of w(v). 

In the case of almost periodic functions the order of averaging can be inter- 
changed and, after another change of variable, we deduce: 


(f) - 2 n A p (+ f _ m fij) A_ q (- f + mfij) W [f - (m + n) fi,]. 
Similarly we obtain for 

k P (r) = E t [8 i (t + r) a (t)] 


(35) 
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the relation: 


r-jfr+jfmfi, r 1 

E t |_e k p (r)J =<!>(+{) A p (- f) 

where 

0 (f) = E r [e-i ft S i (t)], 

By applying (35) and (36) to(18)-(19) and taking (20)-(28) into account we deduce, 
instead of integral equations, the following relations in the frequency domain: 


+ A +1 (+ f - fl x ) A +1 (- f + Dj) Wj (+ f - 2 fij) 

- A +1 (+ f +n,) A_ j (- f -fij) Wj (+ f) 

- A_, (+ f - fij) A +1 *(- f + fij) Wj (+ f) 

+ A_ x (+ f + flj) A_ t (- f - fi,) Wj (+ f + 2 Qj) 


2 7T j S COS 


(- f) -A +1 (- f)], 


(37) 


+ A +2 (+ f - 2 fl,) A +2 (+ f - 2 fij) W 2 (+ f - 4 fij) 

- A +2 (+ f + 2 fij) A_ 2 (- f - 2 fij) W, (+ f) 

- A_ 2 (+ f - 2 fij) A +2 (- f + 2 fij) W 2 (+ f) 

+ A_ 2 (+ f + 2 flj) A_ 2 (- f - 2 fij) W 2 (+ f + 4 fij) 

= T — 1 — - — 0 (f) [A +2 (- f) - A_ 2 (- f )] V 

4 77 j s sin i 1 1 


(38) 


where A p (± f + m flj) is the Fourier coefficient of the term exp [ j(± f + m n x )t] 
in the trigonometrical expansion of a (t). 


19 



If the impulse responses w x (t) and w 2 (t) degenerate into the delta-fiihctions 

W 1 ( T ) = W l 8 (T - Tj), 

w 2 ( T ) = w 2 S (r - r 2 ) (39) 

the integral equations become 

^10 ( r ) “ s w i ( T > T i) cos 1 = 0 . 

20 ( T ) ~ s w 2 ^22 ( T ’ r 2 ) sin i = 0 ( 40 ) 

For the Fourier transforms we have: 


Wi (f) = 


2 77 


and 


W 2 (f) = 



-j f r 

€ 


2 


and the relations (3 7) -(3 8) become: 

Wj {[A +1 (+ f -^ 1 )A +1 (- f + n i )e +2in > T! 

+ A-i (+ f +fi,) A-i (- f 

- [A +1 (+ f + Q t ) A., (- f - Qj) + A_ x (+ f - Q x ) A +1 (- f + fij)]} 


and 


+ j It. 


S J COS 


W„ 


^( f ) [A-i (- f) - A +1 (- f)] , 

{[A +2 (+ f - 2 fi x ) A +2 (- f + 2n i )e +4i ° lT2 
+ A_ 2 (+ f + 2 fij) A_ 2 (- f - 2 n } ) e 4j ^ Tz ] 


(41) 
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- [a + 2 (+ f + 2 flj) A_ 2 (- f - 2 n t ) 

+ A_ 2 (+ f - 2 a,) A + 2 (- f + 2 fi,)]] 


+ j IT 2 

— : : 0 (f) [A +2 (- f) - A_ 2 (- f)] (42) 

sjsini 2 2 ' ' 


These equations show the existence of relations between the Love numbers, time 
lags and Fourier coefficients in the expansion of the tidal effects. 


We deduce the integral equation for w Q (t) by minimizing the quadratic 
functional: 

J' = E |8 fift) - s j_3a o * cos j 

/ +j fl,t\ / -j .0. t\ 

+ a +1 * « J - a_j * ^Wj e J (esc i - 2 sin i) 


" 2 [ a+ 2 * ( w 2 c 21 1 ) ~ a - 2 * ( w 2 e 2jnit )] cos ij | 2 

which represents the mean quadratic error between the observed and computed 


tidal perturbations in Q. After some easy transformations and taking (17*) and 
(26)-(28) into account we can rewrite the last equation in the form: 


J f +00 

k 00 ( T ) w 0 (r) d 

n 


/* +00 +C0 

a oo ( T " w o ( T ) w o W dTd 

•'o -'o 


+ 9 s 2 co s 2 i 


- s' 


n CD 

[(cosec i - 2 sin i) 2 (r, v) (r) Wj (z,) 

j 


+ 4 cos 2 i (r, v ) w 2 (r) w 2 (v)] drdv, 


(43) 
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where we set 

k Q0 (r) = E [S fi(t) a Q (t - r)]. 

Taking the gradient of (43) with respect to w and setting it equal to zero, we 
obtain the integral equation for w Q : 


J fCO 

a oo ( T “ v ) w 0 ( v ) d v = 0 


(44) 


If w {v ) degenerates into a delta-function, 


w o ( T ) = w o 8 C T - T 0 ) 


then (44) takes the form 


k oo ( r ) + 3 w Q s cos i a 00 (r - t q ) = 0, (45) 

which can be used to obtain the Love number w Q . By transforming (44) to the 
frequency domain we have 

Ao (- f) A 0 (+ f) W 0 (f) = -L 0 (f) A, (- f) 

2 TT 

or, taking 

\ (- f) = Aq (+ f) 

into consideration, we obtain a simpler relation 


Aq ( f ) w 0 (f)=2-e (f) 

2 7T 


(44*) 
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where 


0(f) = E t [e"j ft 8 n ( t )] . 

The equation corresponding to (45) is: 

w Q A, (f) = e + jfT °^ (f) ( 45 ') 

and it gives the connection between the Love number w , the time lag and the 
Fourier coefficients in the expansions of the tidal effects in the node of satellite's 
orbit. The relations (40) and (45) or (41), (42) and (45') can be helpful in checking 
the existence of Love numbers from the observations of the long period tidal effects 
in the motion of satellites, 

CONCLUSIONS 

The tidal effects represent a superposition of a large number of periodic 
terms. In the case of nearly equal periods the analytical expansion can easily 
establish the precise form of arguments and the relative significance of amplitudes. 

Long period and resonance terms are especially important. For example, 
the rotation of the lunar orbital plane produces a term of 18 years period in 
tidal perturbations of the ascending node of the satellite's orbit. If these 
effects are not properly taken into account they will contaminate the coefficients 
of the zonal harmonics in the geopotential. 

We do not want to oversimplify our statements. The present theory is not 
a proof, or a disproof, of the existence of Love numbers. However, it must be 
stated, at the same time, that no empirical fit, independently of how successful 
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it is at the moment, can be considered as a complete replacement for theoretical 
thought or a final solution. The long term observations will require a more 
subtle approach, approximately along the lines presented here. We have to 
expect that, because of the rotation of the lunar orbital plane, the tidal effects 
in the motion of some satellites will increase. 

These larger values of the tidal effects and the prolonged observations will 
stimulate a further cooperation between theoretical and numerical directions of 
work in the coming years. 
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Table I. r u (t) 


Coefx 10 4 
of cosine 


V 

^7 

Multiples 
F, r D-jT 

of 

r l T 


Period 
in Days 

2 

( 

0 

1 

1 

1 

1 

1 ) 

16.259 

49 

( 

0 

1 

-1 

1 

1 

1 ) 

83.392 

4 
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1 

-1 

1 

“1 

-1 

1) 

40. 660 

1 243 

< 

0 

2 
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2 

2 

~1 > 

il . 75 8 

1 

( 

0 

2 

— 2 

2 

2 

- i ) 

86. 578 

3 

( 

1 

-2 

C 

— 2 

*® i 

l ) 

20. 51 1 

i 

( 

1 

0 

C 

-2 

0 

l) 

2 3.1 06 

93 

( 

1 

0 

c 

0 

0 

1) 

40. 905 

14 

( 

0 

1 

c 

0 

0 

i) 

09. 809 

89 

( 

0 

2 

c 

0 

2 

i ) 

57. 01 9 


Multiples of 


Coefx 10 4 
of cosine 


V 

'CjT 


Df 

Tj 7 0 1 T 

Period 
in Days 

44 

( 

0 

1 

1 

7 

i ”1) 

11. 738 

2388 

( 

0 

i 

l 

1 

1 

-1 ) 

£ 5. 4£9 

4 

( 

0 

2 

0 

2 

2 

1 ) 

16. 296 

6469 

( 

0 

2 

0 

0 

2 -1 ) 

57. 726 

22 

( 

1 

2 

0 

2 

2 -1 ) 

8. 241 

126765 

( 

0 

0 

0 

0 

0 

1 ) 

84. 428 

5 

( 

1 

0 

0 

— 2 

0 - 1 ) 

51. 046 

24 

( 

1 

0 

0 

0 

0 - 1 ) 

20. 774 

5 

( 

0 

1 

0 

0 

0 -1 ) 

£8. 577 

17 

< 

0 


0 

0 

2 

■ 1 ) 

49. 855 



Table n. 


Multiples of 

Coefx1 ° Period 

of sine ^7 F, 7 D,7 r ,7 ( 1,7 in Days 

“2(01111!) 16.259 

-49 (Ol-lli 1) 83.392 

-4 ( 1 -1 1 -1 -1 1) 40.660 

124 3(0 2 C 2 2-1) 11.758 

1 ( 0 2 -2 2 2 -1) 86.578 

-3(1-2 C ”2 ■■* 2 1) 20.511 

~1 ( 1 0 C -2 0 1) 2 3.1 06 

“53 ( 1 0 C 0 0 1) 40. 905 

-14 (01 C 0 0 1) 09. 809 

-89 ( 0 2 C 0 £ 1) 57.019 


S u( T ) 


Multiples of 


Coefx 10 4 
of .sine 


V 

'CjT F,r 

D,r 


n,7‘ 

Period 
in Days 

44 

( 

0 

1 

1 

1 

1 

- 1 ) 

11. 738 

2388 

( 

0 

1 •= 1 

1 

1 

-1 > 

e5.4fi9 

-4 

( 

0 

2 

0 

2 

2 

1 > 

16. 298 

6489 

( 

0 

2 

0 

0 

2 

- 1 > 

57. 736 

22 

( 

1 

2 0 

2 

2 

> 

8. 241 

126785 

( 

0 

0 

0 

0 

0 

1 > 

84. 428 

5 

( 

1 

0 

0 - 

2 

0 

-1 > 

51 . 046 

24 

( 

1 

0 

0 

0 

0 

-I ) 

20. 774 

5 

( 

0 

1 

0 

0 

0 

- 1 ) 

68. 577 

17 

( 

0 

3 

0 

0 

2 

-1 > 

49.855 



NASA— GSFC 


Table mi 


Multiples of 

Coefx 10* Period 

of sine V ^ l T F 1 T °1 T r i T n i T in Da y s 

-22252 ( 0 2 C 2 2 -2) 10.321 

- 53 168 ( 0 2 C 0 2 ”2) 34.286 

-3(1 2 C 0 2 -2) 15.277 

- A 32 (1 2 C 2 2 -2) 7.508 

-5(2 2 C 2 2-2) 5. SOI 

-6(1 0 C — 2 0-2) 29.1 02 

-3(1 0 C 0 0 -2) 16.672 

-535 ( 0 1 -1 A 1 -2) 42.476 

-1(1 1-1-1 1 -2) 26® 696 

- 152 ( 0 3 C 0 2 -2) 31.345 


S 22 ( r ) 


Multiples of 

Coefx 10* Period 


of sine 


V 

V 

F l^ 

D)T 

r f 


in Days 

-7 

( 

0 

2 

0 

4 

2 

- 2 ) 

6 . C 75 

-6 

( 

0 

2 

2 

2 

2 

- 2 ) 

42 . 745 

17 

( 

1 

-2 

0 

-4 

-2 

2 ) 

7 . 7 S 3 

45 

( 

1 

-2 

0 

-2 

- 2 

2 ) 

16.502 

5960 

( 

0 

0 

0 

0 

0 

2 ) 

42. 214 

66 

( 

1 

0 

0 

0 

0 

2 ) 

79. 348 

-31 

( 

0 

1 

1 

1 

1 

- 2 ) 

1 0 . 3 05 

1 2 

( 

1 

- 1 

1 

-1 

“ 1 

2 ) 

76. 423 

-5 

( 
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1 

0 

0 

2 

- 2 ) 

37.640 



